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In the confined Bethe-Salpeter (BS) methods, the quark propagator is an entire function. Hence, 
the imaginary part of the handbag diagram disappears, leading to a problem of the vanishing parton 
distribution function (PDF). In contrast, the direct calculation of the light-cone (LC) momentum 
distribution does give a non-vanishing result even in the confined BS framework. We consider their 
precise relation and difference, and propose to use the latter as a practical approach to this problem. 
Our formalism is general enough to be applied to various effective models. 

PACS numbers: 11.10.St,13.60.Hb,12.39.Ki,12.40.Yx 



The parton distribution function (PDF) is defined 
through the factorization procedure in the perturbative 
QCD (pQCD) as a low energy-scale quantity which can- 
not be calculated within the framework of pQCD itself 
. It provides us with a unique tool to extract the infor- 
mation on the quark-gluon structure of hadrons directly 
based on experiments. In order to calculate PDF theo- 
retically, the information of the hadronic wave function 
plays the essential role. Hence, it is necessary to resort 
to various low-energy nonperturbative methods, such as 
low energy effective models, the light-front quantization, 
lattice QCD, etc. 

The Bethe-Salpeter (BS) method is one such method, 
which makes it possible to evaluate a number of hadronic 
low-energy observables in a relativistically covariant 
manner 0, 0, 0. The BS equation has been applied 
to mesons which consist of mainly a quark and an an- 
tiquark H H H @. The BS description of baryons 
as the relativistic three quark system has been devel- 
oped in the last decade, utilizing the Faddeev method 
for the quantum three body problem 0, Q, IE @ ■ The 
Faddeev method transforms the three-quark BS equa- 
tion into the relativistic Faddeev equation, which jus- 
tifies the use of the quark-diquark BS equation as its 
approximation|l Ea [H E Hi . 

In the BS method, the effect of the color confinement 
can be taken into account through the entire function na- 
ture of the nonperturbative quark propagator @,0,ELE|- 
Thereby, many of the shortcomings can be improved such 
as the unphysical threshold for decay into colored par- 
ticles, description of higher resonances, etc. However, 
there is actually a serious drawback in calculating PDF. 
Note that, to calculate PDF in the BS method, one usu- 
ally considers the handbag diagram first [Tj, [lfl HH 113 • 
Then, PDF is realized as the Bjorken limit of its imag- 
inary part. In most of the low energy effective models, 
this procedure has been adopted providing reasonable an- 
swers. In the confined BS framework, due to the entire 
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function nature of the propagator, the imaginary part 
of the handbag diagram disappears, leading to a serious 
problem of the vanishing PDF. 

A straightforward solution of this problem may be 
found by considering the intermediate "hadronic loops" , 
which can provide imaginary parts due to the color- 
singlet nature of hadrons. If we adopt this as the so- 
lution, we have to consider the full forward Compton 
amplitude instead of the handbag diagram. Remember 
that, in order to calculate the full Compton amplitude, 
one has to adopt the recently developed "gauge method" 
O fli^ . l20| , which determines the unique set of Feyn- 
man diagrams associated with the particular choice of 
BS interaction kernel. (Other choice of Feynman dia- 
grams lead to the violation of Ward identity.) It follows 
that, in order to obtain the hadronic loop, one has to 
solve a highly complicated BS equation. This is of course 
formidable. Its numerical solution is practically impos- 
sible to be obtained. Moreover, the construction of the 
BS equation itself is extremely non-trivial pit l22j | . The 
transparent relation to the parton picture will be lost as 
well. 

In this way, one has to figure out another practical ap- 
proach to PDF in the confined BS framework. Actually, 
there is another method for PDF by resorting to the di- 
rect calculation of the light-cone (LC) momentum distri- 
bution. Since it is essentially the wave function squared 
with the plus component of the LC momentum of one of 
the partons fixed, it is not expected to vanish unless the 
wave function itself vanishes. This procedure has been 
mainly adopted in the light-front framework, where the 
LC momentum distribution is naturally expressed as the 
equal light-front time correlator. On the other hand, in 
the effective models (including BS method), this proce- 
dure has been scarcely adopted except for a few authors 

mm. 

The latter method is conceptually sound in the 
sense that the Bjorken limit should not be consid- 
ered within the framework of the low energy effective 
theories. The product of the currents should be fac- 
torized at the level of pQCD into the product of the 
high energy scale quantity (the coefficient functions) 
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FIG. 1: The diagram for the matrix element of twist-two 
operators A n . The blob followed by the triple-line denotes 
the BS amplitude for the nucleon, and the single-line denotes 
the propagator of the quark. The cross symbol "x" represents 
the insertion of 7+ (k + ) n ~ 1 /(p + ) n . 



and low scale quantity (the composite operators, i.e., 
the LC momentum distributions). The BS method 
should be used to calculate the matrix element of 
the low energy quantities. Then, the structure func- 
tion F 2 (x,Q 2 ) is expressed for large Q 2 as F 2 (x,Q 2 ) ~ 
x {lu{x, Q 2 ) + \u{x, Q 2 ) + ld(x, Q 2 ) + l~d(x, Q 2 ) + ■ ■ ■ ) 
with the Q 2 -evoluted LC momentum distributions 
u(x,Q 2 ), u(x,Q 2 ), d(x,Q 2 ), d(x,Q 2 ), etc. Fortunately, 
in most of the cases so far, these two methods do not 
result in a serious inconsistency. Especia lly, if the 
propagator satisfies the scaling property 20], both 
methods lead to exactly the same result. This seems to 
be the reason why the latter has been scarcely adopted 
in the effective models. However, the results of these 
two methods do not agree at all in the confined BS 
framework. 

We begin by applying the inverse Mellin transforma- 
tion analytically to the matrix elements of the twist-two 
operators. The aim of doing so is twofold. The first aim 
is to ensure that the direct calculation of the LC mo- 
mentum distribution gives the equivalent result with the 
one in the moment space. As is well-known, pQCD es- 
tablishes the equivalence between PDF in the moment 
space representation and PDF in the x-space represen- 
tation, where the equivalence is provided by the Mellin 
transformation. The second aim is to introduce the quan- 
tities such as the "forward Compton-like amplitude" and 
the 11 handbag-like diagram", which will be used to con- 
sider the precise differences between the two methods for 
PDF. We will argue their physical meaning later. 

We consider the spin averaged matrix element of the 
twist-two operators as follows: 



s=± 



(N(p, s)\ipj+(id + ) n - 1 ip\N(p, sj) 
(p+)" ' 



(1) 



where ip and if) denote the quark fields, and \N(p, s)) de- 
notes the state vector for the nucleon with momentum 
p and helicity s. The light-cone variables are defined as 
p ± = -^7j(p° ±P 3 ), etc. We adopt the covariant normal- 
ization (N(p,s)\N(p',s')) = 2E N (p)(2ir) 3 5 ( - 3) (p-p/)5 ss > 
with En(jj) 



y 



p 2 . 



The actual calculation of 



the matrix elements of these bilinear operators in the 
BS framework amounts roughly to the diagram depicted 




FIG. 2: The contour of the integral in Eq. @. 

in Fig. ^ The cross " x " indicates the point where we 
insert the following vertex as 



7 +(fc+) f 



(2) 



(For more details, see where the matrix element 

of the conserved current is discussed based on the gauge 
method |l8lll9l 20], which may work as explicit examples. 
Ref.0 is also helpful.) 

We apply analytically the inverse Mellin transforma- 
tion to A n . For this purpose, we introduce an analytic 
function T(x) for ieCin the following way: 



1 \ 

K ' 2m ^ x n 

n— 1 



(3) 



Since T(x) is the analogue of the forward Compton am- 
plitude T l j, l/ (x,Q 2 ) in pQCD, we refer to T(x) as "for- 
ward Compton-like amplitude" . We will argue its phys- 
ical meaning later. To proceed, we have to make two 
assumptions on T(x). (1) The power series in Eq. @ 
converges, if |a; > 1. (2) The singularity structure of 
T(x) is the cut along the segment [—1,1]. Due to the 
assumption (1), in the region \x\ > 1, T(x) defines an 
analytic function, which can be analytically continued 
to the inside as much as possible. The most important 
property of the forward Compton-like amplitude T{x) is 
the following identity: 



An — 



dx x 71 - 1 



T(x), 



(4) 



c 



where the contour C is a sufficiently large circle centered 
at the origin as depicted in Fig. [21 This identity follows 
from the residue theorem. Due to the assumption (2), 
we can deform the contour C to the thin contour C as 
depicted in Fig. [21 We are thus lead to the following 
identity: 



dx x n - 1 H{x), 



(5) 



where H{x) is the discontinuity of T(x) along the cut, 



i.e. 



H(x) = disc(T(x) 

= T(x — ie) — T(x + ie). 



(6) 



3 



A 

p .JL 



k 

JL p 



p 




FIG. 3: The handbag-like diagram, which corresponds to the 
forward Compton-like amplitude T(x). The dashed line rep- 
1 7 + 

resents ■ — . The meanings of the blob followed by 

2m xp+ — k+ 

the triple-line and the single-line are the same as in Fig. 



FIG. 4: The diagram which represents H(x). The symbol 
"®" represents the point, where r y + 5(xp + — k + ) is inserted. 
The meanings of the blob followed by the triple-line and the 
single-line are the same as in Fig. Q 



We define the quark and the anti-quark distribution func- 
tions for x (0 < x < 1) as 



q(x) = H(x), q{x) = -H(-x), 
respectively j^. Now, Eq. J^J reads 



^4n. — 



dxx r 



1 (q(x) + (-l) n q(x) 



(7) 



(8) 



From the relation to the moments, this identity shows 
that q(x) and q(x) really work as PDF. We note that 
H(x) is directly expressed in the canonical operator rep- 
resentation in the following way ptij : 



H(x) 



-V / e lxp z (9) 

x(N(p, s) \lp(0h + ^(z-)\N( P , S )), 



where ip(z~) is a shorthand notation for ip(z + — 
0, z~ , z±_ = 0). Indeed, by inserting Eq. © into Eq. JSJl, 
we obtain Eq. (JTJ. 

To see explicitly that q(x) and q(x) in Eq. Q cal- 
culate the LC momentum distribution, we complete the 
summation in Eq. pj analytically. Since the cross sym- 
bol in Fig. n represents the vertex 7+ {k + ) n ~ 1 / (p + ) n , the 
summation in Eq. reduces to the summation of these 
vertices as 



V(x) 



y> 1 7 +(fc+)"- 1 

n— 1 v 7 

1 7 + y> / fc + 

n=0 v ^ 



(10) 



2ni xp 
1 



7" 



2ni xp + — fc+ 

The diagram, which corresponds to the direct calcula- 
tion of T(x), is depicted in Fig. El where the dotted line 
represents V(x). Since this diagram is the analogue of 
the handbag diagram, we refer to it as 11 handbag-like di- 
agram" . In order to calculate H(x), all we have to do is 
to replace V(x) by the following quantity: 

V(x-ie)-V(x + ie) (11) 

= 7^ f 1 

2iri \xp+ — k+ — ie xp+ — k+ + ■ 

= j + 5(xp + — k + ), 



where we used the identity: 5(x) = ^ - . 

Since Eq. I|ll|) counts the "number" of each parton which 
carries the LC momentum xp + , q(x) and q(x) in Eq. Q 
calculate directly the LC momentum distribution. By 
construction, it is clear that the quantity which is equiv- 
alent to the twist-two moments under the Mcllin transfor- 
mation is the PDF obtained as the LC momentum distri- 
bution rather than the PDF obtained from the handbag 
diagram. 

We consider the physical meaning of the handbag- 
like diagram and the precise differences between the two 
methods. We may think of the handbag-like diagram 
as the "proper Bjorken limif of the handbag diagram. 
Indeed, the LC momentum distribution is obtained as 
the imaginary part of this "proper Bjorken limit" of the 
handbag diagram, while, in the original method, the 
imaginary part of the handbag diagram in the straightfor- 
ward Bjorken limit is considered. (We emphasize again 
that these two methods would agree, if the propagator of 
the struck quark could satisfy the scaling property |20j|.) 
We can find the precise difference between the two meth- 
ods in the difference between these two diagrams. In the 
handbag-like diagram, the free propagator is used for the 
struck quark, while, in the handbag diagram, the non- 
perturbative propagator is used. Intuitively, since the 
extremely large momentum is transfered to the struck 
quark, there is no reason to believe in the effective non- 
perturbative propagator, which is constructed to repro- 
duce the low energy properties. Furthermore, since the 
origin of the propagator of the struck quark in our for- 
mulation is the derivatives in the twist-two operators, it 
should not be the non-trivial one. 

Several comments are in order. First, the nontrivial 
propagator should be used except for the struck quark. 
In Ref . ^3 , the first attempt to calculate PDF for pion in 
the confined BS framework is performed, where the free 
propagator is used not only for the struck parton but also 
for the spectator parton. Their calculation is different 
from our formalism in this way. Second, by construc- 
tion, our LC momentum distribution naturally satisfies 
the various sum rules, i.e., the fermion number sum rule, 
the momentum sum rule, and so on. (Sec Eq. JSJ. Its 
moments calculate the matrix elements of various con- 
served currents.) The third comment is concerning the 
calculability of the LC momentum distribution. Usually, 
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to evaluate the LC momentum distribution, one insert 
r ) + 8{xp + — k + ) as the vertex in the diagram like Fig.^ 
In some cases, it is not easy to perform this procedure 
straightforwardly. (For instance, when the integration 
momentum becomes complex due to the Wick rotation, 
the meaning of the delta function with the complex vari- 
able is highly nontrivial.) In such cases, the technique of 
the inverse Mellin transformation can extend the calcula- 
bility of the LC momentum distribution [2j| . As is often 
the case, the calculation in the moment representation is 
much easier than the one in the ^-representation. The 
last comment is concerning the gluon. Once the gluon is 
attempted to be explicitly taken into account, the deriva- 
tive in Eq. QJ is replaced by the covariant derivative. At 
the canonical operator level, the inverse Mellin transfor- 
mation leads to Eq. where the path-ordered integral 
Pexpi J c A^dx^ 1 is inserted. Here, C is the light- like 
straight path connecting and (z + — 0,z~,z± = 0). 
At the BS level, it becomes nontrivial to perform an- 
alytically the Mellin transformation. We note however 
that, it is not the handbag-like diagram with its struck 
quark propagator replaced by the non-perturbative one 
— it would be much more complicated. In this sense, 
the straightforward solution may look attractive, since it 
makes us to avoid this difficulty. However, as we have 
argued, its practical implementation is formidable. An- 
other solution may be provided by A + = gauge. 

To summarize, we have considered the problem of van- 



ishing imaginary part of the handbag diagram in the con- 
fined BS framework. We have seen that a practical ap- 
proach is provided by the direct calculation of the LC 
momentum distribution. In order to consider the pre- 
cise differences between the two methods for PDF, we 
have applied analytically the inverse Mellin transforma- 
tion to the expectation value of the twist-two operators. 
We have introduced the forward Compton-like amplitude 
and the handbag-like diagram, which are the analogues 
of the forward Compton amplitude and the handbag di- 
agram, respectively. In our formalism, the handbag-like 
diagram serves as the "proper Bjorken limit" of the hand- 
bag diagram, and its imaginary part calculates the LC 
momentum distribution, i.e., the PDF. We finally remark 
that our formalism is general enough to be applied to var- 
ious effective models. In realistic calculations, we often 
encounter with the nontrivial quark propagators, which 
cannot satisfy the scaling property. Even if such a non- 
trivial propagators are involved, the evaluation of PDF 
can be done with our method. 
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